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Abstract 

Comprehensive analysis on the light Higgs scenario of the minimal supersymmetric stan- 
dard model is performed in the framework of the non-universal Higgs mass model (NUHM). 
The NUHM is known to be the simplest framework realizing the light Higgs scenario with 
the unifications of gaugino and sfermion masses at the scale of grand unified theory being 
consistent. All important constraints from collider experiments, flavor-changing processes, 
and cosmological observations are considered in order to clarify the allowed region of the 
model parameter space, where the Markov Chain Monte Carlo method is used to obtain the 
region. We show that the region is excluded by current LHC results for the SUSY particle 
searches. We discuss possible extensions to save the LHS. Even with such extensions, the 
measurement of Bs — > fJ-^ f-i^ would enable us to test the LHS in the near future. 
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1 Introduction 



In spite of tremendous successes of the standard model (SM), both hierarchy and dark mat- 
ter (DM) problems require us to consider new physics beyond the SM. The hierarchy problem 
is essentially given rise to by quadratically divergent corrections to the Higgs mass term, so that 
the new physics is expected to appear at around or less than 0(1) TeV scale. On the other hand, 
recent cosmological observations such as the Wilkinson Microwave Anisotropy Probe (WMAP) 
experiment [l][2] have established the existence of dark matter, which cannot be accounted for 
in the framework of the SM. Among several candidates for dark matter, the weakly interacting 
massive particle (WIMP) is one of the most promising candidates |3]. Since the mass of the 
WIMP is predicted to be less than 0(1) TeV scale in order for its correct abundance observed 
today [Tl[2], the WIMP is also likely to be obtained by the new physics at the TeV scale. 

The minimal supersymmetric standard model (MSSM) [4] is thought to be the most promis- 
ing candidate for the new physics. It gives solutions to both hierarchy and dark matter problems 
simultaneously and enables us to unify SM gauge interactions. The direct search of the Higgs 
boson at the LEP experiment severely constrains the parameter space of the MSSM, because 
the model predicts the mass of the lightest neutral Higgs boson h to be the same as or less 
than the Z boson mass at the tree level. In principle, there are two prescriptions to avoid the 
constraint. One is that h is composed mainly of the up- type Higgs. Radiative corrections from 
scalar top quarks then push up the mass of /i [6]. This scenario has frequently been discussed 
so far. In the limit where superpartners are very heavy (although it needs fine-tunings), signals 
of the Higgs boson at collider experiments becomes very similar to that of the SM. 

Another prescription is that h is composed mainly of the down-type Higgs. In this case, there 
are regions to evade the constraint though radiative corrections from scalar bottom quarks do 
not push up the h mass enough. The LEP constraint is evaded because the coupling of the ZZh 
interaction can be suppressed [7HT4]. This is called the light Higgs boson scenario (LHS). The 
LHS can be tested at the LHC experiment because the difference from the SM inevitably exists; 
for example, the masses of other Higgs bosons in the LHS as well as that of h are predicted to be 
around the Z boson mass. The LHS has some attractive features. The LEP anomaly [S], which 
is the excess of the Higgs-like events with the Higgs mass being around 98 GeV reported by 
the LEP collaborations, can be explained. Severe fine-tunings on the supersymmetric (SUSY) 
little hierarchy problem is not required [10^. It is also worth noting that the LHS is consistent 
with a simple framework of grand unified theory (GUT) and seesaw mechanism 15 . The 
cosmological aspects have also studied [TTl[T31[Tni[I7]. Several recent papers have studied the 
expected parameter space in MSSM [18]. 

In the LHS, the charged Higgs boson significantly contributes to the 6 — > 37 process, since 
their mass is very light as that of Z boson. This contribution turns out to be too large, 
because 112^+ > 350 GeV is required for the type II two Higgs doublet model [TH]. This should 
be hence canceled out with SUSY contributions. Because of this reason, the LHS predicts 
that both chargino and stop are light enough. This fact inevitably leads to large Br{Bs — 
/x^/ii"). As a result, the LHS must have large SUSY contribution to 6 — 57, while small 
contribution to Bs — > fJ-^ fJ- - In order to calculate these processes accurately, the LHS should 
be studied comprehensively in the framework of a unification model. In this paper, we analyze 
the LHS in the framework of the non- universal Higgs masses model (NUHM), which is the 
simplest model realizing the LHS in the MSSM. Using the Markov Chain Monte Carlo (MCMC) 
method, we search the parameter region consistent with all important constraints from collider 
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experiments, flavor-changing processes, and cosmological observations. In particular, we impose 
the constraint of the dark matter relic abundance observed by WMAP [1] and the bound by 
direct detection experiments [20]. We find that the LHS in the framework of the NUHM is 
marginally inconsistent with SUSY particle search in the LHC experiments [211 [22] . 

This paper is organized as follows. In the next section, we briefly introduce the LHS in the 
framework of the MSSM, and also explain how the NUHM realizes the LHS. In Section [31 the 
framework of our analysis for the comprehensive study of the LHS is shown in addition to some 
detailed explanations of experimental constrains used in our analysis. All of our physics results 
is presented in Section [U Section [5] is devoted to summary. We also discuss some relaxations 
of experimental constraints and the boundary condition of the NUHM. 

2 Light Higgs boson scenario 

The MSSM predicts two Higgs doublets denoted by iJ„ and Hd- One Higgs doublet i?„ gives 
the masses of up-type quarks, while down-type quarks as well as charged leptons acquire their 
masses by interactions with Hd ■ Five (real) scalar components in Hu and Hd remain as physical 
states; one CP-odd pseudo-scalar Higgs boson A, a pair of charged Higgs bosons ff^, and two 
CP-even Higgs bosons h (lighter) and H (heavier). In the basis of (ReiJ^, ReiJ°) with H^ and 
being neutral components of the two Higgs doublets, the mass matrix is given by 



Here, uia (mz) is the mass of the pseudo-scalar Higgs boson (the Z boson) and {sp) = 
cos/3 (sin /3), where the ratio of two vacuum expectation values (Hu) and (Hd) defines the angle 
/3 as tan /3 = {Hy) / {Hd) ■ Radiative corrections in each component of the mass matrix are 
summarized in A^j , and their detailed expressions are found in Ref . 4 . The mass eigenstate of 
CP-even Higgs bosons h and H are obtained by diagonalizing this matrix. With the use of the 
mixing angle a, the mixing matrix to diagonalize the mass matrix is given by 



The LEP collaborations have given constraints on masses and couplings of the Higgs bosons 
not only in the SM but also in the MSSM [13]. According to their results, the lightest Higgs 
boson h in the MSSM can be lighter than 114.4 GeV , which is the lower bound on the SM 
Higgs boson, when the coupling between h and Z bosons gzzh is suppressed enough. This 
is because the discovery of the Higgs boson at the LEP2 experiment relied on the process 
e+e^ — > Zh through the s-channel exchange of the Z boson. In the MSSM, this coupling is 
given by gzzh = g'zzh sin(/3 — a) with g^'zzh being the coupling between Z bosons and the Higgs 
boson of the SM. In order to obtain the Higgs boson h lighter than 114.4 GeV , a suppressed 
sin(/3 — a) is required. Since tan /3 is expected to be large enough to satisfy experimental 
constraints as shown in the following sections, a large mixing angle a ^ 7r/2 is needed. On the 
other hand, the LEP boimd on the SM Higgs boson is also applied to the heavier Higgs boson H. 
This constraint is complementary to that on h, because the coupling between H and Z bosons 
is given by gzzH = g^zzh cos(/3 — a). Furthermore, the coupling between Z, A, and h is also 
proportional to cos(/3 — a), so that the search of Higgs bosons at the LEP2 experiment using 
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the e^e — > Ah process should be paid attention to. Since the gzAh coupUng originates in a 
derivative couphng, this process does not provide a serious constraint due to P-wave suppression 
when ruh > 90 GeV [l3]. 

Since the constrained MSSM (CMSSM) [23l[23 is too restricted to represent all degrees 
of freedom of the Higgs sector, the LHS is impossible to be realized. Particularly, the soft 
SUSY breaking masses of two Higgs doublet fields mj^ and are fixed by other parameters 
such as universal supersymmetric particle masses mo and to 1/2- As a result, when irih is 
lighter than 114.4 GeV , all other supersymmetric particles become also very light, and such 
a light Higgs boson has already been excluded. The mass spectrum of the Higgs bosons in 
the CMSSM is therefore inevitably the decoupled type, namely, the lightest Higgs boson is 
up- type (ReiJ°, meaning a ~ 0) with the mass of the order of mz, while other Higgs bosons 
are much heavier, of the order of 1 TeV. This fact means that constraints on nij^^ and rn^^ 
should be relaxed, or even be chosen freely from other sparticle masses in order to realize the 
LHS. One interesting possibility is to relax the boundary condition — mj^^ — TOq to be 
= (1 + '^Hu(Hd))"^o GUT scale Mq- Since Higgs multiplets are not necessarily 

the same ones of other matter superfields, these relaxed boundary conditions do not break the 
simple framework of GUTs. The minimal extension to the CMSSM with this relaxation is called 
the NUHM. This model is hence suitable for the search of the parameter space of the LHS as a 
reference model. The NUHM has six model parameters, (toq, toi/2, ^Oj tan/3, /i, to^), where 
(too, TOi/2, ^0) are defined at AIq and others are defined at the electroweak scale. Parameter 
degrees of freedom for fi and toa are translated into mn^ and at Mq ■ This parameterization 
allows us to deal with the masses of two Higgs doublets as free parameters. 



3 Framework of our analysis 
3.1 Markov Chain Monte Carlo 

We use the MCMC method to clarify the allowed region of the LHS. The MCMC method is 
a random sampling algorithm, which gives a series of parameter sets called the Markov chain 
as output [25j . The samples of the chain obey the distribution which is proportional to a 
given distribution function. The distribution function of our interest is a posterior probability 
distribution function of model parameters x under experimental data D. Bayes' theorem tells 
us that the posterior probability distribution function P(x\D) satisfies the following equation, 

_ P{D\x)Pix) 

where P{x) is the prior probability function reflecting our knowledge about model parameters 
X, while P{D\x) represents the likelihood for the distribution function of experimental data D 
at given model parameters x. In our analysis, a linearly flat prior has been used for P{x), where 
P{x)dx gives a constant probability. For experimental data D, we have used constraints shown 
in the upper part of Table [T] Constraints on flDM^i^ and Br(& 57) are assumed to have 
Gaussian distributions. Other constraints are adopted as boundaries of model parameters. The 
last two constraints in Tabled] are applied later to samples generated by the MCMC. 
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Constraints 


References 


Upper bound on gzzh {gzzn) 


See the text 


EE! 


Upper bound on gzAh {gzAu) 


See the text 


EE! 


Relic abundance of dark matter (51dm^'^) 


0.1099 ±0.0062 


n 


Direct detection of dark matter 


See the text 


m 


Anomalous magnetic moment of muon (Aa^) 


9.5 < (Aa^ X IQiO) < 41.5 


m 


Br(5 -> S7) 


(3.52 ±0.25) X 10-'' 


m 


Br(B, ^ 11+ H-) 


< 1.08 X IQ-^ 


m 


Br(B Tv) 


0.68 < r^ssM < 2.76 




The lightest neutralino mass 


> 50.3 GeV 


m 


The lightest chargino mass 


> 103.5 GeV (92.4 GeV) 


m 


Right-handed selectron mass 


> 99.9 GeV (73 GeV) 


m 


Right-handed smuon mass 


> 94.9 GeV (73 GeV) 


m 


Right-handed stau mass 


> 86.6 GeV (73 GeV) 


m 


Sneutrino masses 


> 94 GeV (43 GeV) 


m 


Stop mass 


> 95 GeV (65 GeV) 


m 


Sbottom masses 


> 95 GeV (59 GeV) 


m 


Squark masses (1st and 2nd generations) 


> 379 GeV 


m 


Gluino masses 


> 308 GeV 


32, 


Squark masses (1st and 2nd generations) 


> 1100 GeV 


[111122] 


Gluino masses 


> 750 GeV 





Table 1: Experimental constraints used in our comprehensive analyses. See refs. [2], [33) . and [34] for 
updated constraints on Q,uMh^ , ^Ami ''■iid Br (6 — > S7), respectively. The last two constraints shown in 
the table, which are obtained at the LHC experiment, are used later for the samples generated by the 
MCMC. Values in parentheses are applied only for limited cases. See the text for more details. 



3.2 Constraints used in the MCMC calculation 

Here, we explain some details of constraints used in our MCMC study. In addition to experi- 
mental constraints shown in Table [TJ we also take a theoretical constraint into account in our 
study, which relates to the stability of the vacuum we are living. Constraint from the p parame- 
ter is not imposed because we found that it is always satisfied very well by our MCMC samples 
of the LHS where Higgs bosons have similar masses. 



• Constraints on the ZZh and ZZH coupling constants 

The upper bound on (g z zni / d'zzh' )^ obtained by the search for the e+e~ — ^ ZHi process 
in the LEP2 experiment [13]. This corresponds to the upper bounds on sin^(/3 — a) for Hi — h 
and cos2(/3 - a) for Hi = H in the MSSM. The bound at 95% C.L. is shown in Table 14 
of Ref. [13] as a function of the T-Li mass. Bounds UB;, and UB,- are available for the cases 
of Br(7^i — !> bb) = 100% and Br(7^i tt) — 100%, respectively. For arbitrary values of 
Br(Hi — !> bb) and Br(Hi — J> tt), we use the upper bound UBtot(Hi) defined by 



/ Br(7^i ^ bb)^ Br(Hi — ^^^x 1/2 
V UB^ ^ UB^ 



UB,ot(Hi)^( "^^'"7""^ + "IZ ) ■ (4) 
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In order to satisfy the constraint for the hghtest Higgs boson h, too smah sin(/3 — a) is not 
required in the range 90GeV < mh < 114.4 GeF, and, in fact, sin(/? — a) < 0.5 is enough. On 
the other hand, h can evade the LEP2 constraint even for mh < 90 GeV when sin(/3 — a) is 
smaller than 0.2. The constraint is so tight that such a region is very restricted as mentioned 
in the next section. It is to be noted that another coupling gzzH on cos(/3 — a) inevitably be 
large when the coupling gzzh on sin(/3 — a) is small. Hence, the heavier Higgs boson H tends 
to receive a severe constraint from the LEP2 experiment in the LHS. 

• Constraints on the ZAh and ZAH coupling constants 

Upper bounds on the MSSM suppression factors (cos^(/3 — a) for Hi ^ h and sin^(/3 — a) for 
Hi = H) are obtained also by the search for e+e" — >■ AHi in the LEP2 experiment, which are 
shown in Table 17 of Ref. [13]. Three bounds (we call UBhf,, UBbr, and UBtt) are presented 
there. Bounds UBhb and UBrr are obtained for cases where both of Higgs bosons decay only 
into bb and rf , respectively. The bound UB^t- is given for the case where one of Higgs bosons 
decays into bb and the other does into rf . For arbitrary branching ratios of Hi and A decays, 
we constrain cos'^(/3 — a) and sin^(/? — a) by using UBtot(^i,^) defined by 

/Br(Hi ^Tf)2Br(A^Tf)2 Bv{Hi bbf&riA ^ bbf 
UBtot(rti,A) - I 



Br(Hi ^ Tf)2Br(A ^ bbf + Br(7^i ^ bbfBi{A rff ^ 



(5) 



Since the upper bounds come originally from the production cross section of two scalars with 
the center of mass energy of ~ 200 GeV, we do not have a severe constraint on the coupling due 
to the P-wave suppression when niA + m-Ui 200 GeV. Roughly speaking, the constraint gives 
a lower bound on tua as tua ^ (200 GeV — mh)- 



• Constraint on the dark matter abundance 



The lightest neutralino plays the role of the DM in the NUHM, and its relic abundance ftuMh'^ 
is computed by solving the Boltzmann equation numerically. This abundance is constrained 
by the WMAP experiment [HE]- We have used the five-year result [1] of the experiment 
(Odm/i^ = 0.1099 ± 0.0062) as shown in Tabled which is obtained by the WMAP data alone 
and seems conservative. Recent result is found in Ref. [2] as iluuh'^ = 0.1120 ± 0.0056. The 
relic abundance is calculated with the use of the micrOMEGAs 2.2.CPC code [35] which allows 
us to deal with resonant and coannihilation processes appropriately. 



• Constraints from direct detections of dark matter 



Since all Higgs bosons are predicted to be light in the LHS, the spin- independent scattering 
cross section between the dark matter and a nucleon (a nucleus) is expected to be large. Di- 
rect detection experiments of the dark matter therefore give an important constraint on model 
parameters. The most severe constraint is now given by the XENON 100 experiment [50] • 
We have used the upper bound on the cross section as a function of the dark matter mass in 
our MCMC calculation. This constraint implicitly assumes that the local dark matter density 
around us is 0.3 GeV /cm'^. This constraint is relaxed if we consider a non-smooth distribution 
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of dark matter in our halo. We take into account this uncertainty on the dark matter density 
by relaxing the constraint of XENON 100 by factor 2 [36]. In our work, we set y = for the 
strangeness component in a nucleon [37] . 

• Constraint on — 2 

The experimental result of the anomalous magnetic moment of muon is known to be deviated 
from the SM prediction. The deviation has been reported as Aa^ — (25.5 ± 8.0) x 10~^° [26] 
(See also Ref. [33]), and we use this value with the 2a region, 9.5 < (Aa^ x 10^°) < 41.5, in 
order to constrain model parameters. 

• Constraint on the 6 — > S7 process 

In our MCMC calculation, we have used Br(6 sj) = (3.52 ± 0.25) x 10"'' [27]. Recent ex- 
perimental bound, Br(6 57) = (3.55 ± 0.24 ± 0.09) x lO""*, is found in Ref. [31]. The SM 
prediction to the branching ratio is, on the other hand, Br(6 — )■ 57) = (2.98 ± 0.26) x 10"'' [38] 
which is somewhat smaller than the observed value. Interestingly, the LHS is possible to resolve 
the discrepancy, because the contribution of makes Br(fe — > S7) increased from the SM pre- 
diction [in]. Since all Higgs bosons are as light as the Z boson mass in the LHS, contributes 
to Br(6 — > 57) significantly and this contribution is, in fact, even too large. Actually, in the 
type II two Higgs doublet model, m/f+ > 350 GeV is required [19]. Hence, it should be canceled 
by SUSY contributions [10], for example, by large A-terms. We have computed the branching 
ratio with the SusyBSG code [3^, in which the next-to-next-to- leading order calculation in 
QCD is involved as in the calculation of Ref. [40] . 

• Constraint on the Bg /^^M process 

Diagrams mediated by neutral Higgs bosons contribute to Bs /i^/i^, and their contributions 
are proportional to {ta.n /3)^/m\ [41]. Since neutral Higgs boson masses are of the order of 
mz in the LHS, the branching ratio of this process gives an important constraint on model 
parameters. Actually, the parameter region with tan^ > 10 is almost forbidden as can be seen 
in the next section. We use the current experimental bound Br(Bs n~) < 1.08 x 10~^ 

at 95% C.L. ^28j which is given by combined analysis of results in the LHCb and the CMS 
experiments. 

• Constraint on the B^ t^v process 

The light charged Higgs boson contributes to the B^ t^v process as in the case of Bs — > 
The experimental result on the process is Br(i3* — > T^i^)oxp = (1.63 ±0.39) x 10"'' [27] , 
On the other hand, the SM prediction is Br(B± t±j/)sm = (1-01 ± 0.26) x 10"" [29] . 
somewhat smaller than the experimental value. We define their ratio as r'^'^P = Br(i?^ — >■ 
exp/Br(i3 T*i^)sM- As a result, the deviation between experimental and theoretical 
results is quantitatively given by r'^'^P = 1.62 ± 0.57. Although the uncertainty still be large 
at present, the constraint is potentially important to restrict the LHS, because the additional 
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contribution from increases the deviation 1121 • We have used the formula 



2 



^MSSM ^ Bt{B^ t=^;/)mssm ^ K _ \ tan^ , 



Co 



Br(S± ^ r±i^)sM 1 \'rn]j± J 1 + eotan^ 



ml\n{mljml) , rn? ln(m? /m?) 



(rn? — m? m? ) (ml — m? )(m? — m? , 

V g Ul/\ Ul bn' ^9 bn'^ bn ' 



(7) 



It is worth noting that t-'^SSM always less than one. In our analysis, the experimental result 
with the 2a region, 0.68 < r"''? < 2.76, is applied as the cut for r^SSM_ 

• Constraints from sparticle searches 

In our analysis, we first derive the boundary condition at the GUT scale Mq using tan/3, /i, 
and niA at the mz scale through the renormalization group running which is evaluated by the 
ISAJET 7.75 code [13 ■ Masses of sparticles are then obtained by running back from Mq to mz 
using the renormalization group running again. There are several lower bounds on the masses 
of sparticles as shown briefly in Table [T] More details on the bounds are as follows: 

The lightest chargino : m^± > 92.4 GeV (m^± - m-^^ < 3GeV or < 300 GeV) 

m± > 103.5 GeV (else) 
Charged slepton : m^^ > 73 GeV (^i^ — < 10 GeV) 

mg^ > 99.9 GeV, m^^ > 94.9 GeV, m^^ > 86.6 GeV (else) 
Sncutrino : > 43 GeV {mc, — ra^^ < 10 GeV) 
mp>94GeV (else) 
Stop and sbottom : m^^ > 65 GeV (m^^ — m^^ < 10 GeV) 
m^^ > 59 GeV (m^^ - m^, < 10 GeV) 
m^^ , m^^ > 95 GeV (else). 

Among those, the bound on m ± is particularly important and gives a severe constraint on 
model parameters. This is because it determines the lower bound on toi/2, and then the lower 
bound on the dark matter mass is given through the GUT relation on gaugino masses. On the 
other hand, the bounds on squark masses also give important constraints on model parame- 
ters because large |Ao| and/or fi tend to be favored in the LHS. Constraints on squark mass 
{rriq > 379 GeV) and gluino mass {rrig > 308 GeV) used in the MCMC study are pre-selections. 
We use further cuts, rriq > 1100 GeV and rrig > 750 GeV, for generated MCMC samples later 
in order to satisfy recent results at the ATLAS and CMS [55] experiments. 



• Constraint from the vacuum stability 



Too large ^-terms gives rise to charge and color breaking (CCB) minima. We therefore take 
the constraint \At\'^ < 3(r7i?^ + m?^ + m'j^^ + /i^) into account in our MCMC calculation. The 
constraint guarantees that the vacuum we are living is the absolute minimum even if there 
are CCB minima [48 . Since rnj^^ + /i^ ^ which are much smaller than stop masses to- 
and ni^^, the constraint gives an upper bound on \At\ in terms of the stop masses. The above 
constraint can be weaker when we also consider the case that our vacuum is metastable whose 
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(GeV) 

Figure 1: Distribution of 4.0 x 10^ samples in the (m?i, sin a)-plane obtained by the MCMC calculation, 
where 2.2 x 10^ samples with |sina| > 0.9 are adopted as the LHS samples in following analyses. 

lifetime is much longer than the age of the universe. In such a case, the constraint is relaxed 
to be + 3//^ < 7.5(m~^ + w?^) [JS]- We adopt, however, the first one to constrain model 
parameters, because it is less ambiguous than the weaker constraint. It is also worth noting 
that the problem of CCB minima occurs when is much larger than mo, and such a param- 
eter region is not favored in the view point of the naturalness. The extension of the allowed 
parameter region to larger Aq is, hence, not so attractive. 

3.3 LHS samples 

With the use of the constraints discussed in Sec. 13.21 we have generated 4.0 x 10^ samples in the 
range rrih < 114.4 GeV through the MCMC algorithm. The distribution of the samples on the 
{mh, sin a)-plane is shown in Fig.[T] Two distinct regions can be seen in the figure. In the region 
around m/j = 114.4 GeV, the lightest Higgs boson h is composed almost of the up-type Higgs 
boson Rei?° because of small |sinQ;|. Thus, the physics property in the region is essentially 
the same as the one frequently discussed in comprehensive analyses of the CMSSM [53]. On 
the other hand, in the region with large | sina|, h consists dominantly of the down-type Higgs 
boson Rei/°, which is characteristic of the LHS as discussed in Seed] We therefore use only 
2.2 X 10^ samples satisfying | sina| > 0.9 in the next section. 

4 Physics results 

There are six parameters in the NUHM. However, only two constraints on JIdm^^ and Br(fe — > 
sj) are used as Gaussian distributions in our MCMC calculation. Central values of those 
constraints are reproduced simultaneously not only at a point but in a wide parameter space. 
When we see number distributions of samples on two parameter space, the region of large 
number of samples does not always means better agreement with experimental data than other 
regions; the region of large number of samples may mean just the large spread of the acceptable 
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region in the projected four parameter space. Therefore, it is better to use the hkehhood in order 
to see the experimentahy acceptable region of the parameter space. The hkehhood quantifies 
the agreement with experiments, and we define it as the product of two Gaussian distribution 
functions for ftuuh'^ and Br(6 — > 57). The normahzation is defined so that the set of the central 
values for them in Table [T] makes the likelihood unity. The larger value of the likelihood means 
the better agreement with experiments. The maximum value of the likelihood is calculated in 
each cell in two parameter space, and we show the distributions in this section. In the region 
where the maximum of the likelihood is close to unity in the two parameter space, experimental 
constraints can be satisfied very well by appropriate choices of other four parameters. 

4.1 Parameter region 

First, results for the parameter distribution in the NUHM are shown in Fig. [2] as the correlation 
of mi/2- The relations, (a) toi/2 - ?7t-o, (b) mi/2-fJ-, and (c) mi/2-^0 are shown, respectively. 
Other two parameter, ttia and tan/3, are distributed in restricted regions, toa ~ 94-107GeV 
and tan/3 ~ 7-10, respectively. The yellow (bright) region is of a large value of likelihood 
maximum in a cell. The contribution to 5 — >■ 57 from charged Higgs should be cancelled out 
by SUSY contribution, and hence allowed parameters have upper bounds. Allowed parameters 
in Fig. [2ja) are divided into two regions by dark matter annihilation processes. One around 
= 140 GeV is pseudo-Higgs mediation region, and the other is coannihilation region. In 
the former region, the relic abundance of the LSP is governed by the pseudo-Higgs mediated 
process, XiXi A ^ bb. In comparison between these two regions, the likelihood is somewhat 
better for the coannihilation region. 

As seen in Fig. [2lb), smaller value oi fj, < 400 GeV are entirely rejected. This is due to the 
direct detection constraint by XENON 100. Since the LHS has the small Higgs masses and a 
large tan /?, direct detection rate which is given by mediation of the lightest Higgs (down-type 
one) is predicted to be large. For a small fi, the rate becomes too large because of additional 
enhancement due to a significant mixing between bino and Higgsino. 

In the coannihilation region, the appropriate relic abundance of the LSP is given by a tuning 
of the masses of the LSP and stau, i.e. a tuning of mi/2 and mg. Since the relic abundance 
in the region is independent of the value of /i, a broad range of is allowed. In the pseudo- 
Higgs mediation region, the relic abundance depends on the value of the term as well as the 
mass difference between the LSP and pseudo-scalar Higgs boson mass because the coupling of 
the LSP to Higgs bosons depends on ^. Meanwhile the relic abundance is independent of toq- 
Large /i is favored in order to satisfy the constraint by direct detection experiments. On the 
other hand, as concerns little hierarchy, too large value of fi is not favorable because light Higgs 
masses require the cancellation between very large values of fi and tti^^^ . 

In Ref. [IH], the authors pointed out that there is a region where the relic abundance is 
determined by the processes XiXi ~^ hA, HA,W^ , Zh processes for lower value of ~ 
300 GeV. For this small /i, since the mixing between bino and Higgsino is significant, these 
processes can be dominant for the control of the relic abundance of the LSP. However, due to the 
smallness of /i, this region is excluded by the recent constraint of direct detection experiments. 

In Fig.[2Kc), the distribution of Aq is shown. The relic abundance of the LSP is independent 
of ^0, and hence the allowed region is almost determined by the 5 — >■ S7 process and the 
condition to avoid dangerous vacua. The large negative Aq is favored by the b ^ sj process 
in order to cancel a large contribution while too large |^o| leads to dangerous charge and 
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Figure 2: The distribution of input parameter correlation. The relations, (a) mi/2-'Tio, (b) mi/2-A'i 
and (c) mi/2-Ao are shown, respectively. The yellow (bright) region is of a large value of likelihood 
maximum in a cell. i n 




color breaking minimum. As a result, the allowed range of has upper and lower bounds. In 
the coannihilation and the pseudo-Higgs mediation regions, \Ao\ becomes large for a large mi/ 2 
and mo, respectively. This trend compensates heavy SUSY particle masses for the contribution 
to & — S7 process. In the pseudo-Higgs mediation region, SUSY particles have larger masses 
than those in the coannihilation region. Therefore, a larger \Aq\ is required in the pseudo-Higgs 
mediation region than that in the the coannihilation region for rni/2 — 140 GeV. 

In addition, the Bg fJ-^ fJ'~ process gives complementary condition to the & — 57 process 
and the relic abundance. The branching ratio of Bg A*"*"/^" is proportional to tan^ (3/m\. 
Hence, this process is drastically enhanced in the LHS due to a small itia- As a result, large 
tan /3 > 10 is ruled out by this constraint. On the other hand, the lower bound on tan/3 comes 
from the realization of the LHS 8 . The lower bound tan/? > 7 results in Br(£?s ji^ > 
7 X 10^^ (See Fig. [3]). Constraint on this process also excludes the region where lighter stop 
mass is very light as 200 -300 GeV due to a large \Aq\. A light stop mass is favored for h s^, 
and however this inevitably enhances Br(i?s ^'^fj.~). 

Before the discussion about the LHC constraint on SUSY particles, we mention the distri- 
butions of Higgs boson masses. The SM-like Higgs search at the LHC constrains the mass of the 
heavier neutral Higgs in the LHS because it is composed mainly of the up- type Higgs. Although 
stop masses are not large in LHS, the SM-like Higgs mass is lifted up because of following two 
contributions; one is from the large A-term which is required to be consistent with the 6 — > S7 
constraint and the other is from the off-diagonal component of the neutral Higgs mass matrix. 
As a result, the mass is found to be distributed in the range of 114 GeV < mn < 122 GeV in 
our MCMC samples. Since the LHS requires a large |sina| (namely, ttia ^ mz), the allowed 
region of is nearly the same range; 94 GeV < < 107 GeV. Because of = m\ + 
at the tree level, is restricted as 122 GeV < mjj+ < 133 GeV. The ranges of mH+ and 

tan /3 are not constrained by present results of the charged Higgs search by collider experiments. 

Finally, let us consider LHC constraints on SUSY particles. In the analysis, we use the 
NUHM as a benchmark model in which the constraint on boundary conditions of mn^ and 
niHd are relaxed as compared with CMSSM. But the constraints on the parameter space of 
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(mo, m.1/2) would be exactly similar to the CMSSM case in the LHS region shown above. In 
these models, differences arise only from mn^ and mu^ values which can be rewritten by 
and /i at the weak scale while the masses of gluino and first generation of squark are nearly 
the same. Additionally, the LSP mass is also nearly the same in the LHS region shown above 
because the Higgsino component of the LSP is small. Then, the production cross section of 
SUSY particles and collider signal of these events would be no change and the constraint on the 
CMSSM parameter space pTl[22] excludes the whole LHS region shown above. 

5 Summary and Discussions 

We have analyzed comprehensively the LHS of the MSSM in the framework of the NUHM, 
which is the simplest one realizing the LHS in GUTs. We have taken into account, all important 
constraints from collider experiments, flavor-changing processes, and cosmological observations. 
To search allowed parameter sets, we have used the MCMC method. We have shown that there 
are two parameter regions consistent with these constraints except for LHC results. One is the 
coannihilation region and the other is the pseudo-Higgs mediation region. To reconcile the 
contribution to the b ^ s"f process with experimental bounds, all SUSY particles remains light, 
as < 1000 GeV. Hence, we have found that the current LHC results have already excluded the 
whole LHS region in the NUHM. 

Finally, we discuss that some relaxation of our restriction or some extensions of our simplest 
model could have allowed LHS regions in the rest of this paper. First, we consider the relaxation 
of experimental constrains on JIdm and muon g — 2. Other constraints are considered to be 
robust or have not reduced sample points signiflcantly. Second, the restriction of the NUHM 
is relaxed. That is, the boundary condition at the GUT scale is relaxed, or the gauge group of 
the GUT is broken at an intermediate scale above the EW scale. 

The relaxation of the restriction of the relic abundance of the LSP is expected to bring us 
the large number of allowed parameter sets. It is favorable that the thermal relic abundance 
of the LSP is consistent with the WMAP observed dark matter abundance. However, if this is 
not the case, the model is not excluded unless the relic abundance is larger than the WMAP 
abundance. The relic abundance can significantly be reduced for a smaller /i, since the Higgsino- 
bino mixing is enhanced. The reduced abundance of the LSP relaxes the constraint from direct 
detection experiments. Then, a small /i which enhances the direct detection cross section 
is allowed in this relaxation condition. In space (or irig-mq space), the extended 

allowed region of /i does not give new sample points at the neighborhood of the coannihilation 
regions where the relic abundance is independent of fi. However, a new allowed region with 
H < 400 GeV appears for mi/2 ^ 300 GeV. In this region, the LSP is significantly annihilated 
by XiXi ~^ HA, W^H^ , Zh processes and larger values of mo and mi/2 are allowed. Then, 
the squark and gluino masses can exceed the LHC constraints. 

Next, we mention the relaxation of the muon g — 2 constraint. In pseudo-Higgs mediation 
region, rriq is extended to rriq ~ 1.2 TeV by removing the g ~2 constraint. However, this relaxed 
region is still rejected by the LHC results because the gluino mass is light as 400 GeV. On the 
other hand, the coannihilation region is extended to mi/2 — 450 GeV. For this upper limit, 
the gluino mass is about 1000 GeV. This remains as the region excluded by the LHC results. 
Hence, the relaxation of the g — 2 constraint does not bring us new allowed parameter sets. 

The relaxation of NUHM conditions changes the spectrum of the MSSM particles. If the 
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gauge group of the GUT is broken into some gauge group (for example, SU(3)^ x SU(2)^ x 
SU(2)^ X U(l)g_^ of the left-right model) at an intermediate scale above the EW scale, the 
mass relation of the MSSM particles can be modified. In the modified spectrum of the MSSM 
particle masses, first and second generation squarks are required to be heavy because of the 
LHC constraints while stop should be light in order to have a sizable SUSY contribution to 
b — >■ 57. This might be realized through the renormalization running with a large |^o|- Another 
possibility might be that the third generation squarks are very distinguished from first and 
second generations as a boundary condition. Of course, the MSSM without any boundary 
conditions at the GUT scale would have enough parameter space to realize the LHS allowed by 
current experimental data. Detailed studies of each possibilities are beyond the scope of this 
paper. 

Even in new parameter sets obtained by such relaxations, Higgs bosons should be light in 
the LHS and they give large contributions to 6 — > 57 and Bg — > /-i^^i^. Gaugino and stop should 
also be light enough to cancel a large contribution of to 6 — > 57. Furthermore, the charged 
Higgs boson production is enhanced because the coupling, gw^HTh^ depends on cos(/3 — a) [Sj 
(For the case with CP-violating scalar sector, see also Ref. [50] and references therein). As a 
result, the LHS can be judged in the near future by the LHC SUSY particle search, the LHCb 
search of Bs — > M^/^^ process and the search for H^h production at the LHC experiments. 

Note added: Updated constraint on Bs l-i-^l-J-^ at the CMS just appeared as Br{Bs — > 
/i+^~) < 7.7 X 10^^ at 95% C.L. [51 during finalization of this article. A small region in 
Fig. [3] satisfies the bound, and the LHS would be tested in near future without assuming mass 
spectrum of SUSY particles. 
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